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I.  INTRODUCTION 


This  report  discusses  an  adaptive  algorithm  for  optimizing 
the  gain  of  an  antenna  array  on  an  incoming  signal.  The  algorithm 
is  based  on  a  steepest-ascent  maximization  of  the  array  output 
power,  subject  to  a  constraint  on  the  array  weights. 


This  algorithm  is  under  investigation l for  use  in  conjunction 


ioruft 


with  a  power  inversion  adaptive  array -E4-5-J,  a  modified  version  of 
the  LMS  adaptive  array  A  power  inversion  adaptive  array  can 

provide  significant  protection  from  interference  in  a  spread  spectrum 
communication  system,  but  it  does  not  provide  any  form  of  beam  tracking 
on  the  desired  signal.  Beam  tracking  can  be  obtained,  however,  by 
first  arraying  antenna  elements  in  pairs  with  power  inversion  feed^ 
back,  and  then  combining  the  element-pair  outputs  with  the  gain 


optimizing  algorithm  described  here.  -4n-this  report,  wp  discuss 
only  the  gain  optimizing  algorithm. The  combined  system  will  be  V 
the  subject  of  a  future  report. 


n 


izU**-5 


II.  A  GAIN  OPTIMIZING  ALGORITHM 

Consider  an  array  of  N  antenna  elements  as  shown  in  Fig.  1. 

For  simplicity,  the  elements  are  assumed  to  be  isotropic  and  non¬ 
interacting.  The  signal  from  each  element,  Yi ( t) ,  is  passed  through 
a  processor  Pi  that  generates  K  outputs  labeled  X-i (t)  on  Fig.  1. 

Each  output  X-j (t)  is  multiplied  by  a  real  weighting  coefficient  Wf 
and  then  is  summed  to  produce  the  array  output  S(t).  The  processor 
Pi  may  be  either  a  quadrature  hybrid  with  two  outputs  as  shown  in 
Fig,  2,  or  a  tapped  delay-line  with  K  outputs,  as  shown  in  Fig.  3, 

A  quadrature  hybrid  processor  (Fig.  2)  provides  a  simple  magnitude 
and  phase  adjustment  of  the  signal  Y-j ( t ;  and  is  the  appropriate  form 
of  processing  when  the  signals  are  narrowband.  A  tapped  delay-line 
processor  (Fig.  3)  provides  a  frequency  dependent  transfer  function 
behind  each  element  and  is  appropriate  for  wider  bandwidth  signals* 
The  gain  optimizing  algorithm  to  be  developed  in  this  report  may  be 
used  with  either  type  of  processing,  and  an  example  of  each  will  be 
given  below. 


In  this  report  we  develop  an  iterative  algorithm  for  adjusting 
the  weights  Wi  such  that  the  array  gain  is  maximized  on  an  incoming 
signal.  The  algorithm  is  based  on  a  steepest-ascent  maximization  of 
the  array  output  power,  subject  to  a  constraint  to  prevent  the  weights 
from  going  to  infinity.  The  algorithm  discussed  here  is  in  the  spirit 


♦Reference  [8]  contains  a  quantitative  comparison  of  the  bandwidth 
performance  of  quadrature  hyprid  and  tapped  delay-line  processors 
for  a  two-element  array. 
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Fig,  3.  Array  with  tapped  delay-line  processing. 


of  the  LMS  algorithm  developed  by  Widrow,  et  al .  [6],  and  discussed 
in  many  reports  on  adaptive  arrays  [7,9].  The  general  properties  of 
constrained  gradient  techniques  have  been  discussed  by  Eveleigh  [10] 
and  by  Gill  and  Murray  [11]. 

Let  us  assume  that  each  of  the  signals  Xi(t)  is  sampled  periodic¬ 
ally  in  time,  with  an  interval  between  samples  equal  to  At.  Let  X-i ( j ) 
denote  the  value  of  X^(t)  at  the  jth  sampling  instant  and  also  let 

r*,(«  l 


X(j)  - 


"1  / 
x,(j) 


[xNK(j)J 

be  a  column  vector  with  components  X-j(j).  Also,  let  us  assume  that  at 
the  jth  sampling  instant,  the  weights  in  the  array  have  the  value  Wi(j), 
which  may  be  represented  by  the  column  vector  W(j); 


We  wish  to  develop  an  iterative  algorithm  by  which  these  weights 
can  be  adjusted  at  each  sample  to  maximize  the  array  gain  on  an  incom¬ 
ing  signal . 

Intuitively  speaking,  maximum  array  gain  on  an  incoming  signal 
will  result  in  the  greatest  output  signal  power  from  the  array.  Since 
the  array  output  signal  is 


(3) 


NK 


S(t)  =  l  W.X.(t) 
i=l  1  1 


the  average  output  power  is 


NK  NK 


(4)  S‘(t)  =  l  l  WfVO.MX.U) 

i=l  1=1  1  3  1  J 


where  the  overbar  denotes  the  time  average.  We  would  like  to  choose 
the  weights  so  S&(t)  is  maximum.  However,  there  is  clearly  no  upper 
limit  to  S2(t)  if  the  W-j  are  allowed  to  take  on  arbitrary  values,  so 
it  is  necessary  to  maximize  S^(t)  subject  to  some  constraint  on  the 
values  of  the  weights.  Several  types  of  constraints  appear  to  be 
possible,  but  in  this  report  we  assume  the  weights  must  satisfy  the 
constraint  equation 


(5) 


NK 

l 

1=1 


Wi2  -  1 


I.e.,  the  weights  are  constrained  to  lie  on  the  surface  of  a  hypersphere 
of  unit  radius. 

To  optimize  the  gain  of  the  array,  we  adopt  the  following 
iterative  algorithm 


(6) 


W(j+1)  =  W(j)  +  kv(j) 


where 

W(j)  =  weight  vector  at  the  jth  iteration 

W(j+1)  =  weight  vector  at  the  (j+l)st  iteration 

k  =  scalar  constant  controlling  the  rate  of  adaptation 

and  v(j)  =  a  vector  correction  term  chosen  to  move  the  weights 

toward  the  maximum  gain  setting. 

The  correction  term  (j)  is  obtained  by  taking  the  gradient 
of  S*(t)  and  retaining  only  its  component  parallel  to  the  hypersphere. 
Since  the  gradient  of  S'^ ( t )  is  given  by 


we  find  from  Eq.  (4)  that 
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This  gradient  vector  can  he  written  as  the  sum  of  two  components  -- 
one  perpendicular  to  the  weight  hypersurface  and  one  tangent  to  the 
hypersurface.  That  is, 


(9)  vw(S2)  =VW((  (S2)  +  (S2) 

where  vw  (S  )  and  vw  (S  )  are  tangent  and  perpendicular,  respect¬ 
ively,  to  the  surface 'Of  the  hypersphere,  as  shown  in  Fig.  4.  To 
obtain  vw^  (57),  we  let  denote  a  unit  vector  normal  to  the  surface 

of  the  hypersphere,  (n  is  also  shown  in  Fig.  4.)  n  will  be  given  by 

w  w 

no)  fiw  = 

where 

(ii)  R  =  J  Wi2  +  w22  +  .  •  .  wNK2 

2 

vw  (S  )  is  then  given  by  (T  denotes  the  transpose) 


W1 

W 

w„ 


w, 
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(12) 


Vw  (s2)  =  [C  7w{s2)]  K 


W^S  +  W2X^S  + 


WNKXNKS 


A 


W. 


W, 


NK 


SURFACE  OF 
THE  HYPERSPHERE 


<S2) 


Fig.  4.  Perpendicular  and  tangential  components  of  v  (S£). 
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By  subtracting  ^  (S2)  from  vw(S2)  (see  Eq.  (9)),  we  obtain 

X^S  S2!^ 

R2 

(13)  vw  (S2)  =  :  -  s\ 

n  — rk 


In  view  of  the  constraint  equation, 
~  NK  9 

(14)  R2  =  I  W  2  =  1 
i=l  1 

2 

we  may  drop  the  R  term  to  obtain 


\  ($2)  ■ 


S(X1  -  w^) 
S(X2  -  W2S) 


S^XNK  '  WNKS* 


Since  we  would  like  to  control  the  weights  in  real  time,  it  is  impractical 
to  compute  the  time  average  of  the  quantities  S(Xi  -  W-jS)  indicated  by 
the  overbar  in  Eq.  (15).  It  would  be  possible  to  compute  the  average 
of  this  quantity  over  a  finite  time  period  and  consider  this  as  an 
estimate  of  the  infinite  time  average.  However,  as  is  done  in  the  LMS 
algorithm  (6),  we  will  adopt  the  simplest  estimate  of  all,  namely,  we 
ignore  the  averaging  completely  and  just  use 


"n  - 


S(X,  -  W,S) 
*  S(X2  -  WZS) 


l  S<XNK  -  WNKS>J 

The  iterative  algorithm  given  in  Eq.  (6)  thus  becomes 


(17)  W.(j+1)  =  W.(j)  +  kS(j)  [X.(j)  -  W.(j)S] 

Let  us  first  consider  what  happens  with  this  algorithm  if  the 
weights  are  off  the  constraint  surface  for  some  reason.  We  would  like 
the  algorithm  to  be  such  that  if  the  weights  drift  off  the  surface, 
they  will  automatically  return  to  the  surface.  With  this  algorithm, 
it  turns  out  they  will  do  this.  To  see  why,  note  that  the  second 
term  In  Eq.  (13)  represents  the  perpendicular  component  of  the 
gradient.  In  simplifying  Eq.  (13),  we  obtained  Eq.  (15)  by  dropping 
the  term 

,  NK  9 
RZ  =  l  W  2  , 
i=l  1 

which  is  unity  if  the  weight  constraint  is  satisfied.  Note,  however, 
that  if  the  weights  are  off  the  hypersphere  for  some  reason,  the 
equation 


8 


will  not  hold.  Hence,  the  second  vector  in  Eq.  (15)  will  be  parallel 
to  the  perpendicular  component  of  the  gradient,  but  will  have  the 
wrong  magnitude.  If,  for  example, 

l  w/  >  1, 

the  vector 


has  a  greater  magnitude  than  vw^  (S'j.  Then  in  Eq.  (15),  in  the 
second  term,  we  subtract  a  vector  greater  than  the  perpendicular 
component  of  vw(S^),  Note  that  since  the  perpendicular  component  of 
vw(S*)  always  points  outward  away  from  the  origin,  (we  can  always 
increase  the  output  power  S?  from  the  array  simply  by  increasing  all 
weight  coefficients  in  proportion),  the  second  term  subtracted  in 
Eq,  (15)  overcorrects  for  the  perpendicular  component  and  the  resulting 
total  vector  in  Eq.  (15)  points  back  toward  the  hypersphere;  i.e., 
it  has  the  correct  parallel  component  but  has  a  perpendicular  compon¬ 
ent  pointing  inward.  This  inward  perpendicular  component  will  cause 
the  weights  to  move  back  onto  the  constraint  surface. 

In  a  similar  way,  if 

I  w,2  <>. 

the  second  term  in  Eq.  (15)  will  have  a  smaller  magnitude  than  the 
perpendicular  component  of  vw(S^).  The  net  vector  in  Eq.  (15)  will 
not  have  its  perpendicular  component  completely  cancelled— there  will 
be  a  residual  radial  component,  which  will  move  the  weights  back  out 
toward  the  surface. 

Thus,  we  see  that  dropping  the  term 

l  w,2  ■  ' 

in  going  from  Eq.  (13)  to  Eq.  (15),  which  was  done  to  simplify  the 
algorithm,  is  important  because  it  makes  the  weights  always  tend 
toward  the  constraint  surface. 
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Finally,  we  discuss  the  stability  of  the  iterative  algorithm 
in  Eq,  (17).  Because  the  signals  are  sampled  at  time  increments  At 
and  a  correction  is  applied  to  the  weights  at  each  sample  instant, 
the  system  controlling  the  weights  may  be  viewed  as  a  sampled  data 
control  system.  Hence,  we  may  expect  that  if  the  loop  gain  constant 
k  is  too  large,  the  system  will  become  unstable. 

To  obtain  an  idea  of  the  suitable  range  of  values  for  k,  let  us 
suppose  that  in  the  steady  state,  the  fth  array  weight  will  have  the 
value  W0^.  (Note  that  the  final  steady-state  solution  is  not  unique  — 
there  are  many  sets  of  array  weights  that  will  maximize  the  array 
gain  toward  a  given  signal.  Woi'  represents  one  possible  set  of  steady- 
state  weights.)  Suppose  at  the  jth  iteration,  W.  (j)  differs  from 
WQi  by  AW.(j):  1 

(19)  W.(j)=Woi+  AW.(j) 


Using  this  equation  in  Eq.  (17)  yields 


(20)  WQi  +  AW.(j+l)  =  WQ.  +  AW.(j)  +  kS( j )[X. ( j)  -  W0.S(j)-AW.(j)S(j)] 

Since  Wpi  is  a  steady-state  solution  to  the  weight  iteration  equation, 

W01*  will  have  the  property  that  the  average  value  of  the  quantity 


S(j)  [X.(j)  -  w0.s(j)] 

is  zero;  i.e.,  in  the  steady-state,  the  correction  term  in  Eq.  (17) 
will  average  to  zero  over  many  samples,  so  both  Wi(j)  and  Wi(j+1)  assume 
the  same  value,  W0i.  Cancelling  the  W0i  term  and  dropping  the  term 
S(J)  [X.(j)  -  WoiS(j)]  from  Eq.  (2)  leaves 

(21)  AWi(j+l)  =  AW.(j)  -  kS(j)AW.(j)S(j) 

-  AW.(j)  [1  -  kS2(j)] 


In  order  for  the  algorithm  to  converge,  the  average  value  of  the 
factor  1  -  kS2(j)  must  be  such  that 


AW.U+l) 

AW.(j) 


(22) 


|1  »  kS2(j)|  <  1 


so  that  at  successive  iterations  the  difference  between  Wj(j)  and 
Mol  becomes  smaller.  To  satisfy  Eq.  (22),  k  must  lie  in  the  range 

(23)  0  <  K  <  — 1 — 

2S2(j) 

Furthermore,  since 

(24)  S(j)  =  XTW  =  WTX  , 

we  have 

(25)  S2(j)  =  WTXXTW  . 


Note  that  the  constraint  equation 

Z  wi2  - 1 

implies  that  the  vector  W  is  a  unit  vector.  Hence,  the  maximum  possible 
value  of  sZ(t)  will  equal  the  maximum  eigenvalue  of  the  matrix 

(26)  *  =  XXT 


(This  case  will  occur  if  the  unit  vector  W  lies  along  the  principal 
axis  of  $  associated  with  this  eigenvalue.)  Thus,  if  a  denotes 
the  maximum  eigenvalue  of  $,  the  feedback  algorithm  wilnBe  stable 
as  long  as  k  lies  in  the  range 

(27)  0  <  k  <  — L_  , 


We  can  obtain  a  more  readily  usable  upper  bound  for  k  by  noting  that 

NK  — 5— 

(28)  atox  *  Trace  $  =  £  X^(t)  , 

so  the  algorithm  will  be  stable  if  k  is  restricted  to  the  narrower 
range 


I 


tv 


U.:*. 


1 

(29)  0  <  k  <  - 

NK  -5 - 

2  I  Xf(t) 

i=l  1 

Note  that 


NK 


x^t) 


is  proportional  to  the  total  power  incident  on  the  array. 

Examples 

Now  we  give  two  examples  of  the  use  of  this  algorithm.  The 
first  example  is  a  two-element  array  with  quadrature  hybrid  processing 
and  a  narrowband  signal .  The  second  example  is  a  two-element  array 
with  tapped  delay-line  processing  and  a  wide  bandwidth  signal. 


EXAMPLE  1: 

Consider  a  two-element  array  as  shown  in  Fig.  5.  A  signal  is 
assumed  to  propagate  Into  the  array  from  an  angle  relative  to 
broadside.  (Ue  assume  the  antenna  elements  to  be  isotropic  and  non¬ 
interacting.)  We  will  assume  also  that  there  is  no  noise.  As  a 
result,  the  received  signals  In  the  elements  are  given  by 


(30)  Y.j  (t)  ■  a  cos  [«Qt] 


(31)  Y2(t)  =  a  cos  [«Qt  -  ♦q] 


where  a  Is  the  amplitude  of  the  signal,  w0  is  the  carrier  frequency, 
and  ♦  is  the  interelement  phase  shift  due  to  the  propagation  delay: 


02)  *0  *  ~  sin  e 


(L  Is  the  element  spacing  and  A0  Is  the  free-space  wavelength.)  The 
signal  Y-j (t)  Is  arbitrarily  chosen  to  have  zero  electrical  phase 
angle. 


Fig.  5.  A  two-element  array  with  quadrature  hybrids. 


The  signal  from  each  element  is  passed  through  a  quadrature 
hybrid;  the  signals  out  of  the  hybrids  are  denoted  by  X-j ( t) .  X](t) 

and  X?(t)  are  the  in-phase  and  quadrature  components,  respectively, 
of  Y-j(t),  and  X?(t)  and  X4U)  are  the  in-phase  and  quadrature  com¬ 
ponents  of  Y2(t).  Thus,  we  have 


(33)  X^ ( t)  =  a  cos  («0t) 

(34)  X2(t)  =  a  sin  (u)Qt) 

(35)  X3( t)  =  a  cos  (u>0t  -  4>q) 


(36) 


XA(t)  =  a  sin  (u  t  - 


A  Fortran  computer  program  was  written  (see  Appendix  I)  to 
simulate  the  behavior  of  the  iterative  algorithm  in  Eq,  (17).  In 
the  examples  to  be  shown  below,  we  have  arbitrarily  chosen 

(37)  a  =  1 


(38)  L  =  aq/2  (half-wavelength  spacing), 
and  for  the  initial  value  of  the  weight  vector. 


(Note  that  this  satisfies  the  constraint  equation.) 

Some  results  of  this  simulation  ar  displayed  in  Figs.  6  to  13. 

In  Figs.  6  and  7,  the  desired  signal  arrives  from  broadside 
(e  =  0°).  Figure  6  shows  the  transients  that  result  in  the  four 
weights,  and  Fig.  7  shows  the  final  array  pattern  after  the  weight 
transients  have  ended.  It  may  be  seen  that  the  beam  maximum  points 
in  the  proper  direction,  and  it  can  be  shown  that  the  final  array 
pattern  has  the  maximum  possible  gain  in  the  direction  of  the  desired 
signal.*  Figures  8  and  9  show  similar  results  for  the  case  when  the 
desired  signal  arrives  from  9  *  30°.  Figures  10  and  11  show  e  *  60°, 
and  Figs,  12  and  13  show  e  =  90°.  In  all  cases,  the  final  weights 
in  the  array  yield  a  maximum  possible  gain  in  the  signal  direction, 
and  it  can  be  seen  from  the  patterns  how  the  beam  is  steered  toward 
the  signal. 


*For  a  broadside  signal,  maximum  gain  occurs  when  W,  =  W,  and 

if  _  i.i  10 
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140. 
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Fig.  6.  Weight  transients  0=0°  (quadrature  hybrids) 
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EXAMPLE  2 

We  again  consider  a  two-element  array  but  with  a  two-section 
tapped  delay  line  processor  behind  each  element,  as  shown  in  Fig,  14 
Each  delay  line  section  is  a  quarter  wavelength  long  at  the  center 
frequency  of  the  signal. 


x^ttL 


(T)  "« 


"•  (l)  w» 


s  (t) 


Fig.  14.  A  two-element  array  with  tapped  delay-lines. 


An  amplitude  modulated  signal  is  assumed  to  propagate  into  the 
array  and  to  produce  element  signals 


Y-j  (t)  =  a[l  +  cos  u)mt]  cos  [«  t] 


(41) 


Y2(t)  =  a[l  +  cos  <um(t>T)]  cos  [«  (t  -  t)] 


where  a  is  the  amplitude,  <%,  is  the  modulation  frequency,  id0  is  the 
carrier  frequency,  and  t  is  the  propagation  time  delay  between  elements, 
given  by 

(42)  t  =  ”  sine 

Each  delay  line  section  (between  taps)  is  a  quarter  wavelength  long 
at  frequency  u>0  and  hence  causes  a  time  delay  T  equal  to 

<«>  T  -  s£- 

o 

Thus 

(44)  X^t)  =  Y]  (t) 

(45)  X2(t)  =  Y^t  -  T) 

(46)  X3(t)  -  YT(t  -  2T) 

(47)  X4(t)  *  Y2(t) 

(48)  X5(t)  =  Y2(t  -  T) 

and 

(49)  X6(t)  =  Y2(t  -  2T) 

The  feedback  algorithm  in  Eq.  (17)  has  been  simulated  in  this 
problem  (see  Appendix  II)  with  the  following  parameter  values 

(50)  a  =  1 

<51>  *  =r 

(52)  L  =  ^ 
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and  the  arrival  angle  e  *  0°,  30°,  60°,  and  90°.  The  initial 
weight  vector  was  chosen  to  be 


1 

0 


0 

0 


The  weight  transients  and  the  resulting  array  patterns  (computed  at 
frequency  u>0)  are  shown  in  Figs.  15  through  22. 

Figures  15  and  16  show  the  weight  transients  and  the  final 
array  pattern  (at  frequency  <d0)  when  0  =  0°.  Figures  17  and  18  show 
similar  results  for  0  =  30°,  Figs.  19  and  20  show  0  =60°,  and  Figs. 
21  and  22  show  0  =  90°.  In  all  cases,  the  resulting  weight  settings 
maximize  the  pattern  response  toward  the  signal. 


CONCLUSIONS 

This  report  has  discussed  an  iterative  algorithm  for  adjusting 
the  weights  in  an  adaptive  array  to  maximize  the  array  gain  on  an 
incoming  signal.  The  algorithm  is  based  on  a  steepest-ascent  maximiza¬ 
tion  of  the  array  output  power  subject  to  the  constraint  that  the 
sum  of  the  squares  of  the  array  weights  is  constant.  It  was  shown 
that  the  algorithm  prevents  the  weights  from  drifting  away  from  the 
constraint  surface  and  also  that  the  algorithm  is  stable  for  a 
suitable  range  of  the  feedback  gain  constant. 

Two  examples  showing  that  the  algorithm  does  optimize  gain  were 
presented,  one  with  a  CW  signal 'and  guadrature  hybrid  processing 
behind  the  elements,  and  the  other  with  a  modulated  signal  and  tapped 
delay-line  processing  behind  the  elements. 
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Fig.  18,  Final  pattern,  0*30°  (tapped  delay-lines), 


Fig,  22.  Final  pattern,  e  *  90°  (tapped  delay-line) 
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